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^ . The chiral condensate in QCD at zero temperature does not depend on the quark chemical 

' potential (up to one third the nucleon mass), whereas the spectral density of the Dirac 

(N ■ 

operator shows a strong dependence on the chemical potential. The cancellations which make 
this possible also occur on the microscopic scale, where they can be investigated by means 
of a random matrix model. We show that they can be understood in terms of orthogonality 
^ • properties of orthogonal polynomials. In the strong non-Hermiticity limit they are related 

pH ' to integrability properties of the spectral density. As a by-product we find exact analytical 

expressions for the partially quenched chiral condensate in the microscopic domain at nonzero 
, chemical potential. 
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I. INTRODUCTION 

One of the simplest questions one can ask regarding the phase diagram of QCD is: What is 
the dependence of the chiral condensate on the quark chemical potential at zero temperature? 
This simple question has an equally simple answer: At low temperature, the vacuum of QCD 
is dominated by the pions and since they have zero baryon charge, the free energy remains /i- 
independent until /U surpasses a third of the nucleon mass. Despite its simplicity, even this answer 
is exceedingly hard to verify by direct lattice QCD computations. To see why, let us consider the 
relationship between (the magnitude of) the chiral condensate 

S(m) = -—log Z(m) = -( Tr- \ (1) 

and the spectral density p of the Euclidean Dirac operator, P + ij,^q. At nonzero chemical potential, 
the Dirac operator is non-Hermitian and the spectral density spreads out from the imaginary axis. 
Since the chiral condensate is related to this /^-dependent spectral density by 

1 ^,2. 1 



= T7 / — \ p{z,z*,m;fi) (2) 

V Jr z + m 
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it is far from obvious that the chiral condensate obtained this way is independent of /i (for < mjv /3 
and zero temperature). 

The quark mass enters the support of the Dirac spectrum when /i > , and by an 

electrostatic analogy ^ one could conclude that the chiral condensate becomes /i-dependent and 
always vanishes in the chiral limit (m — > 0) when the chemical potential is non-zero. While this 
conclusion is correct for quenched QCD, it fails completely in the unquenched case where, as argued 
above, the unquenched chiral condensate is independent of 

In the unquenched case, the generally accepted picture for two decades was that chiral symmetry 
breaking requires that the Dirac spectrum should accumulate to a non-zero density on the imaginary 
axis (like the Banks-Casher relation \a\ at /_f = 0) despite the fact that the anti-Hermiticity of the 



Dirac operator is explicitly broken by the chemical potential (see [3[ for a discussion of this point). 
Indeed, intricate cancellations due to dynamical fermions set in when /i becomes greater than 171^^/2. 
The eigenvalues, however, do not accumulate on the imaginary axis. Direct input P, 6| from the 
microscopic domain of QCD shows that the discontinuity in the chiral condensate arises due to 
strongly oscillating terms in the eigenvalue density Q]. The samemechanism is also responsible for 



Because of the oscillations, the 



chiral symmetry breaking in one dimensional QCD at nonzero ^ 
unquenched eigenvalue density is not real and positive. This fact is a direct a consequence of the 
sign problem: since the measure of the Euclidean QCD partition function includes a complex valued 
fermion determinant the expectation value of a real and positive function, such as Ylik "^^(-^ ~ ^k), 
needs not be real and positive. 

The microscopic domain Q] is the region where the eigenvalues of the Dirac operator are in the 
domain 

\z\ < i-^, (3) 

while the volume is much larger than Aq^j-,. The quark masses may or may not be in this region. 
In this region the QCD partition function can be expressed in terms of microscopic scaling variables 

z = zVT,, rh = mVT, and fi = fiFj^VV, (4) 

which stay fixed in the thermodynamic limit. This scaling limit is known as the microscopic 
limit. The oscillations in the spectral density with a period proportional to the inverse volume are 
resolved at this scale. The cancellations of these oscillating contributions with an amplitude that 
grows exponentially large with the volume results in a ;U-independent chiral condensate. They give 
us a direct insight into problems faced by lattice QCD at nonzero chemical potential. 
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FIG. 1: Phase diagram of the Dirac spectrum in the complex z plane for one dynamical flavor. The support 
of the quenched spectrum is between the two vertical black lines. Unquenching introduces the oscillations 
within the ellipses which intersect the a;-axis -aX x — ±m and x ~ ±(|F^/i^/I] — m). The labels N, tt. K 
refer to the normal, pion and kaon condensed phases of the generating functional of the eigenvalue density. 
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In the /i-independence of the chiral condensate was established using an asymptotic approx- 
imation to the exact expression for the microscopic spectral density. The argument was carried 
through by complex contour integrations. In the present paper we extend this result to the complete 
mass dependence of the chiral condensate in the microscopic domain. 

The extent of the oscillating region in the complex eigenvalue plane can be understood from 
a mean field treatment of the chiral Lagrangian for the generating functionals of the eigenvalue 
density (see jlo|). The picture that emerges is shown in Fig. |T] where we present results for one dy- 
namical flavor. Inside the green ellipses the eigenvalue density is strongly oscillating. Outside these 
regions the oscillating terms are exponentially suppressed and the eigenvalue density approaches 
the quenched spectral density. The spectral density can be evaluated by introducing valence quark 
masses z and z* in addition to the dynamical quark with mass m. The first two masses can be 
thought of as the up and down quark whereas the mass m can be interpreted as a strange quark 
mass. Varying the chemical potential then corresponds to varying the isospin chemical potential 
and the strangeness chemical potential explaining the appearance of regions we identify as the pion 
condensed phase (denoted by vr), the kaon condensed phase (denoted by K) and the normal phase 
(denoted by N). In the latter phase the chiral condensate is nonzero but the Bose condensates are 



vanishing. For more details we refer to [10|]. 

The fact that the structure of the eigenvalue density is determined by the phases of the gener- 
ating functional shows that this structure persists beyond the microscopic domain. In this paper, 
however, we will focus on the microscopic limit. In this limit the QCD partition function and the 
spectral density of the Dirac operator are uniquely determined by the global symmetries of QCD 
and hence can be derived from chiral random matrix theories with the same symmetries [U|. In 



random matrix theory, the microscopic limit can be taken by identifying the size of the matrices 
with the volume of space-time and defining a scaling limit as in A random matrix model for the 
microscopic domain of QCD at nonzero chemical potential amenable to an analysis by orthogonal 
polynomials was introduced in j^, and its properties were further analyzed in [g] and Q]- 

The purpose of the present paper is two-fold. First we will show that within the random matrix 
model [5|] the delicate cancellations follow from orthogonality properties of orthogonal polynomi- 
als. The results are exact, independent of the size of the matrices. It allows us to establish the 
^-independence of the full microscopic chiral condensate as obtained from integrating the micro- 
scopic spectral density according to Eq. ([2]). We also present a direct computation using complex 
contour integrations which establishes the relation between the exact /i-dependent microscopic 
spectral density and the /^-independent chiral condensate. As a by-product we will obtain the 
exact expressions for the partially quenched chiral condensate at non-zero chemical potential as 
well as the bosonic partition functions in an arbitrary fixed topological sector. 

In the upcoming section the random matrix model is introduced along with the complex orthog- 
onal polynomials and their relation, through the Cauchy transform and kernels, with the eigenvalue 
density. In section IIIII we then use the orthogonal polynomials to establish the /x-independence of 
the chiral condensate. A strategy to derive this result without the use of orthogonal polynomials is 
sketched in section HVl but it could only be worked out in detail in the strong nonhermiticity limit 



(see Appendix C ) . In section |V] we turn to the bosonic partition function and the quenched and 
partially quenched condensates are discussed in section IVII . In the appendices we provide alter- 
nate derivations of the results of section IIIII Appendix A reproduces the microscopic condensate 



starting from the microscopic density so that only universal quantities are used in the calculation. 
Appendix B uses a method involving complex contour integration while in Appendix C we first 



rewrite the eigenvalue density total derivative. 
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II. THE RANDOM MATRIX MODEL 

We consider the microscopic limit of QCD at nonzero chemical potential. In this limit where 
(S is the magnitude of the chiral condensate in the chiral limit and i<V is the tree level pion decay 
constant) 

mSy and ^l^FlV (5) 

are kept fixed for y — > co, the mass and chemical potential dependence of the QCD partition 
function is given by a random matrix model (see for example the review [l^). For quark flavors 
with mass m and n pairs of regular and conjugate quarks with masses x and y*, respectively, this 
partition function is defined by 

xdet"(P(^) + X )det"(pt(^) + y*)^ (6) 

where the non-Hermitian Dirac operator is given by 

(0 i^/a + u,^-i/a\ 

Here ^ and ^' are complex (iV + z^) x matrices both distributed according to a Gaussian weight 
function 

wg{X) = exp(-iVTrX^X) . (8) 

The parameters a and 7 are scale factors used to map the random matrix model onto the chiral 
Lagrangian, which are given below. The number of additional rows as compared to columns 
gives rise to v zero modes of P, and u is therefore referred to as the topological index. Inverse 
determinants are interpreted as bosonic quarks and will be denoted by negative values of iVj or n. 
In cases where Nf or n are zero they will be left off along with the corresponding masses. 
In the random matrix model, the microscopic limit is given by the limit iV — > 00 where 

m = 2Nam and fi^ = 2N-f^fi^ (9) 

are kept fixed as — > 00. The identification with the QCD partition function is made according 
to (see the discussion in Q]) 

m = 2Nam mT.V, (10) 
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This then determines the scale factors a = YW I2N and 7^ = F^V/2N . From here on we will drop 
the factors of a and 7 except when explicitly needed. 

Contrary to the Hermitian random matrix ensembles, it is quite nontrivial to express the par- 
tition function ([6]) as an integral over the joint probability distribution of the eigenvalues of 'D(fi). 
Remarkably, it was shown in [s] that an analytical form could be obtained with result given by 

r ^ 

Z^^'^''{m,x,y*;f,)r^m''''f{xy*r / lldhkV^f'^{{z,},{zn,m,x,y*; (11) 

•^'^ k=i 

where the integration extends over the complex plane and the joint probability distribution reads 

1 ^ 

V''f^^{{z,},{zn,m,x,y*;fi) = — \Ar,{{zf})\' l[w{z,, zl; n){m' - zlff {x' - zlTiy*' - z^'T ■ 

(12) 



k=l 



The Vandermonde determinant is defined as 



N 



i>j=i 

and the weight function includes a modified Bessel function. 



(13) 



w{zk,zl;fi) = \zk\ 



iV(l-/.^) 

2^2 ,^.iyexp^ ^ 



(14) 



The modified Bessel function is obtained as a result of the integration over the angular degrees of 
freedom. In the microscopic limit the same functional dependence is obtained starting from a chiral 
Lagrangian for the bosonic phase quenched partition function [3]. This is a strong argument for 
the universality of this factor. 

The eigenvalue representation makes it possible to evaluate integrals over eigenvalues by means 



of the method of comp 



ex orthogonal polynomials 



16 



13, 



Il8l ]. Analytical expressions for 

the eigenvalue density [5| , eigenvalue correlation functions [sl and partition functions [3, [3] have 
been obtained using this method. 



A. Orthogonal polynomials and their Cauchy transform 

In this section we introduce the orthogonal polynomials and relate them to the quenched and 
unquenched eigenvalue density. The complex orthogonal polynomials corresponding to the weight 
function (jl4p are given in terms of the complex Laguerre polynomials by 

p.(-;/^) = (V^)'^!^^(-Y^)- (15) 
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They satisfy the orthogonality relations 

(fzw{z,z*;fi)pk{z;fi)pi{z;n)* = 6ki r"^ , (16) 



with the norm r'^ given by 



Since the orthogonal polynomials pk are related to the Laguerre polynomials they also satisfy an 
orthogonality relation on the imaginary axis, 



dx Wi{x) pk{ix] fi) pi{ix; fi) = 6ki s^, (18) 



with 



We will make use of this relation in section IIIIB 11 Since 

n._, = PJ^zPl^ (20) 

is a polynomial of order Z — 1 in z^, as a direct consequence of the orthogonality relations, we can 
establish the identity 

d^zw{z,z ■,ii)pi[z ) 5 5 — = 0. (21) 

C z^ — 

The Cauchy transform of the orthogonal polynomials is defined as 

hk{m;fi)= I d'^z^r-^ — -w{z, z* ; n)pl{z; jd), (22) 
z^ - 

where we recall that C indicates that the integration extends over the complex plane. 

The partition function for one fermion can be expressed in terms orthogonal polynomials as 

— = zpN[z;n), (23) 



and the partition function for one bosonic flavor is given by a Cauchy transform [y, llSl . Il9l ] 



(24) 



Note that at finite N both partition functions depend on the chemical potential. This dependence 
can be removed [2^ from the fermionic partition function (j23p by scaling the factors a and 7 in ([7]), 
but we will not bother here since the ^-dependence will factorize in the microscopic limit which is 
our main concern. The bosonic partition function has a nontrivial /i-dependence that can not be 
removed through the constants a and 7 



IS 



21[ . Below we will drop fx from the argument of p and 



h to make the notation less clumsy. 
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B. Kernels and the spectral densities 

Prom the orthogonal polynomials and their Cauchy transforms we can construct two different 
kernels (see for example [a]) for later use 

nj f ^ ^ Pk{x)hk{y) . 
T-CN[x,y) = 2^ . (26) 

k=o 

We will also use the auxiliary kernel 

AfN{x,y)= ^ +nN{x,y), (27) 

y X 

which enters in the computation of the partially quenched chiral condensate in section IVI[ 

The /i-dependent spectral density is related to the kernel /C. In particular, the quenched spectral 
density is given by 

P%{z,z*;fi) =2w{z,z*;fi)ICN{z,z*), (28) 



and the unquenched Nf = 1 spectral density can be written as 

(29) 



Pn^ ^{z,z*,m;n) =2w{z,z*;i2) 



)Cn{z,z*) - ^^j^\ }CN{m,z*) 
PN{m) 



Note that the first term is equal to the quenched spectral density. The second term is responsible 
for the strong oscillations of the real and imaginary parts of the unquenched eigenvalue density 
(see figure [2]). There is also an extra contribution zv5^(z) to the eigenvalue density which arises due 
to the exact zero modes of the Dirac operator which has not been included in the spectral density. 

III. THE //-INDEPENDENCE OF THE MICROSCOPIC CONDENSATE 

In this section we will show how a /i-independent chiral condensate in the microscopic limit 
can arise even though the eigenvalue density of Dirac operator which enters ([2]) has a strong fi- 
dependence. Inserting the finite N one flavor partition function (123p in ([T]) we immediately find 
that the chiral condensate is given by 

yS(m) = ^?^^#/^ + ^. (30) 
PN{m) m 

Making use of orthogonality relations for the polynomials we will show that the above expression 
also follows from a direct integration of the spectral density (|29p . The microscopic result is then 
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FIG. 2: The eigenvalue density of the random matrix model ([6]) with N = 20, ^ = 0.8, m = 0.6 and v = 0. 
The quenched density (top) is real and positive. The real part of the unquenched density (bottom) shows 
oscillations. Note that the oscillating region starts out at z = ±m and extends outward. The peak of the 
oscillations has been cut; the maximum amplitude is an order of magnitude larger than the scale displayed. 

obtained by taking the microscopic limit of the finite N result, and as we shall see explicitly 
below, the ^-dependence drops out. The topological term z^/m follows trivially from the 6 function 
contribution due to the zero modes. 

Alternatively one can start from the microscopic limit of the spectral density and show by 
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explicit evaluation of the integrals that the corresponding chiral condensate does not depend on 



the chemical potential. For completeness, this derivation is given in Appendix A The techniques 
used here will also be useful when calculating the partially quenched chiral condensate in section 

ED 

A. Chiral condensate at finite N 

In this subsection we will calculate the one flavor chiral condensate by integrating over the 
eigenvalue density (j29p using a variant of ([2]) obtained from exploiting the evenness of p in z 

VT.(m) = / (fz — 7^ — 7jpZ^~^(z,z*,m;^i). (31) 
Jq — 

The spectral density for Nf = 1 given in Eq. (j29p can be written out as 

A:=0 



By subtraction and addition of pfc(m) and pN{fn) to pk{z) and pn{z), respectively, we obtain 

(33) 



Pn^ ^{z,z*,m;fi) = 2w{z,z*;p) 



Pk{z*)\pk{z) -Pk{m)] _ pk{z*)pk{m) \pn{z) - PNim)] 



ru ^ Th pM(m) 

k=0 " k=0 ^Jv V / 



The first sum in the square brackets is well behaved in the thermodynamic limit whereas the second 
sum, with the exception of z ~ m, shows oscillations with an amplitude that grows exponentially 
with the volume. The chiral condensate is given by the integral 



yS(m) = 2m / (fz'w{z,z*;fi) 



c 



sr^ Pk{z*)[pk{z) -pkim)] _ ^ Pk{z*)pk{'m) [pn{z) -pAf(m)] 
^ Tkim? - z"^) ^ Tk pN{m){'m? - z^) 



(34) 



The contribution of the first sum in (j34p integrates to zero due to orthogonality. To obtain a result 
that is well behaved in the thermodynamic limit, the exponentially large terms in the second sum 
have to be canceled which indeed happens because the contributions to the chiral condensate are 
localized on z = m. The reason is that, as a consequence of the orthogonality of the p^, the sum 

, . .""^ Pk{z*)pk{m) _ 
w{z,z -p) y (35) 

k=0 " 

is a reproducing kernel (or equivalent to the delta function d^{z — m)) in the space of polynomials 
in z of order less than N (denoted by n„(z)). We thus have 

f <fz w{z, z*-p) y Pfc(^*)Pfc("^) n„(^) = iin{m) for n<N. (36) 
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Using the definition of derivative we finahy find the following result for the chiral condensate, 

2m PN{z)-pN{'m) dpN{m)/dm 
= — -— hm 2 2 = r~\ — • ^^'J 

Clearly this is consistent with expressing the chiral condensate as the derivative with respect to 
m of the logarithm of the partition function (j'23|) . Notice that since the spectral density (|29p does 
not contain the trivial contribution of the zero modes, u5'^{z), the condensate calculated this way 
does not have the u/m term. A second remark is that at finite A^, the polynomial p]\j{m) depends 
on the chemical potential. As will be shown in the next subsection, this dependence drops out of 
the chiral condensate in the microscopic limit. 

B. Microscopic limit of the finite N results 

In this section we discuss the microscopic limit of the argument presented in the previous section. 
First, let us show explicitly that, in the microscopic limit, we obtain the expected /i-independent 
chiral condensate. The microscopic limit of the orthogonal polynomials is given by 

(«,)-j^(2Ara,i)^^£:!^ with t = ./|. (38) 



y/rk 7/^\/7r V A^' 

Because of this factorization, the overall ^-dependence of the orthogonal polynomials drops out of 
equation (j37p in the microscopic limit resulting in a /x-independent chiral condensate 

^ lim H!V2iVa) ^ dUra)/dm _ v_ 
^ ' N->oo 2Na Iy{rh) m ^ ' 

The term i//m is canceled by the contribution of the zero modes. In the above argument we have 
manipulated the non-universal finite A^ expressions into a form which are suitable to take the 



microscopic limit. In Appendix A we will show that it is also possible to take the microscopic limit 
from the start and deal with the universal expressions throughout the argument. 

In fact, a slightly stronger result is valid. In the microscopic limit, orthogonal polynomials of 
the same order but with different arguments have the same /^-dependence (see (I38p ). Therefore, 
the ratio 

PN{m){m'^ - z^) ^ PNym) 

does not depend on the chemical potential in the microscopic limit, and for the same reason, each 
term in the sum does not depend on /U. If we define Yij^ as 

2mpk{m) 

^k{m) = -——Ck{m), (41) 

PN[m) 
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so that 

E(m) = J2 ^k{m), (42) 
fc=0 

then also the are independent of /i in the microscopic hmit. Exphcit results for the microscopic 
limit of 5]fc(m) will be worked out in the next subsection. 

1. Microscopic limit ofT,k 

In this subsection we derive the microscopic limit of Sfc(m) defined in (|4ip . This quantity, 
denoted by Ti{m, t), also enters in derivation of the chiral condensate by complex contour integration 
(see Appendix CD - 



The coefficient (m) defined in ([10|) can be calculated by using the orthogonality of the pk on 
the positive imaginary axis 

Ckijn) = / dxwi{x) 2- 2 • y^^l 

Jo si + m-^ 

Now looking at the first part of the r.h.s. of the previous equation we define 

Ei(,n) = 2,„?^M f, «M_^ (44) 



In the microscopic limit we have 



V ^ /" 2Ntdt with t = ^/kjN and u;i(x) = (45) 

JO 



fc=o 

so that 



T}{m,t) = 2mtrh~^ Iy(rht) I dx^ :r^Jv{xt)=2mtIy{rht)Ky{rht). (46) 

Jo a;^ + 

Likewise the second part of Ck{in) contributes to the chiral condensate as 

vii/ A o ^'fc^"^) ^ I ^Pk{ix)pN{ix) 

^fc("^) = 2?7T' — 7 — ri7 / dxwi{x) ^ ^ — . (47) 

In the microscopic limit this simplifies to 

S"(m,t) = 2mt:^^ /""(ix ., ^ . , J^(xt) J^(£). (48) 
I;y(m) Jo ^ + 

The integral is known 2^ resulting in 



S"(m, t) = 2mt^^^Y^Ky{m)Iy{mt). (49) 
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We conclude that the microscopic hmit of ^^(tti) is independent of the chemical potential. 

Reminding ourselves of the identities 

/■I 1 

Jo ^ 

[ dttK^{mt)Iuimt) = l[K^{m)I^{m) + K^+i{m)I^^i{^^^ (50) 
Jo 2 

we obtain after using the Wronskian identity 

S(m)= [\t[t\m,t) -±'\m,t)] = - ^, (51) 

Jo I,^{m) m 

which is the correct result after adding the contribution of the zero modes. 

We have by now established the /i-independence of the microscopic chiral condensate as ex- 
pressed by an integral over the eigenvalue density. As discussed next, the properties of the spectral 
density that lead to the ^u-independence of the chiral condensate can be exposed further by writing 
the the spectral density as a derivative of the partially quenched chiral condensate. 

IV. THE SPECTRAL DENSITY AS A TOTAL DERIVATIVE 

In section IIIII we have seen that the chiral condensate can be written as an integral over a 
reproducing kernel in the space of complex orthogonal polynomials up to order N — 1: 

S(m) = -2m [ (fzKN{z\m) ^^-}^''} . (52) 
Jc PN{m) 

In the thermodynamic limit the kernel can be written as 

lim KN{z*,m) ~ 2m6^(z^ - n?) = -d^* ^ ^ ^ . (53) 
N-*oo IT z^ — 

Using a partial integration this suggest that the condensate can be written as 



2m 

C z^ — 



S(m)= / dPz ^_ ^ d,*F{z,z*,m;ii), (54) 



or that the spectral density can be written as a total derivative with respect to z*. Such a repre- 
sentation of the eigenvalue density in the complex plane is well known [3, Q] ; the function F is 
known as the partially quenched condensate 

F(...-.»;rt4(Tr^^J-^). (55) 



It is the chiral condensate evaluated at a complex mass z. A partial integration of (I54p singles out 
z = m and consistently reproduces the chiral condensate S(m) = F{z = m, z* = m, m; 
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If we can express the eigenvalue density as a total derivative p = dz*F{z, z* ,m; p) and show 
that F at z = z* = m is independent of fx we have an alternative way to show how a strongly /i- 
dependent density can result in a chiral condensate that does not depend on the chemical potential. 



This alternative proof will be worked out in Appendix C in the strong non-Hermiticity limit for 
arbitrary v. While the exact expressions for the microscopic eigenvalue density are known, the 
partially quenched condensates at nonzero chemical potential have not been obtained previously 
and will be derived in section IVII 



V. CAUCHY TRANSFORM AND THE BOSONIC PARTITION FUNCTION 

In this section we evaluate the Cauchy transform of the orthogonal polynomials for finite N. 
As a new result we obtain the microscopic limit of the bosonic partition function with arbitrary 
topological charge. For zero topological charge the result agrees with the bosonic partition function 

The result 



obtained in 



1^ and thus proves a conjectured cancellation in the derivation of 
of this section can also be used to obtain an alternate expression for the microscopic limit of 
7i]y{mi,m) at finite N. 

The Cauchy transform is defined by (see Eq. (j22p ) 



hk{m) = [ cfzw{z,z*;fi) ^ ^ _y, Pk{z*) 
Jc 

\2v+2 



f ^^K^(6|^|2)e-(-^+-*^)p,(^*) (56) 
z^ -m^ 



ic 

with a = N(l — fi'^)/4:fi'^ and b = N{1 + p?)/2^'^. The integrals can be evaluated by writing 



_j,2 _ ^2 _ j^2 -,2 _ j^2 



+ .2 ^2 • (57) 



The first term has a pole while the second is analytic. The angular integration for the first term is 
now easy since the integrand only depends on z* (or \z\) and not z. Simply expanding l/{z'^ — m?) 
in a geometric series in z^ jm^ we obtain 



\ Sz ^ ^J^,(6|z|2)e-'^("^'+^*')pfc(z*) = / drr2'^+3^,(6r2)e-'^(™'+^'/"^')pfc(rVm) . 

For the second piece, since it is analytic, we assume it can be expanded in the form 

oo 

e-''''"Y.'^k{m)pk{z) . (58) 



Z2 



fc=0 

Then by substituting this into ([56|) and using orthogonality we get 



> e e 



1^1"^^ ^2_^2 ^.(&kP)e-"^* Pk{z*) = rldk{m), (59) 
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where is the normalization factor defined in (jl7p . To evaluate dk{m) we first assume that the 
series (|58p converges for all z and m. Then this is also valid along the imaginary axis z = ix where 
we can use the orthogonality relations of the orthogonal polynomials Pkiix) defined in (jlSp . The 
coefficients are thus given by 

-a{x'^+m^) 



4(m) 



1 



dxwi{x)pk{ix) 



k Jx=o 



Using that 



dt e-*(-'+-') 



-a{x'^+nn?) _ ^ 

x"^ + 



and the formula 



this can be reduced to 
dkiin) = ■ 



(l-/.2)2/4^2 



A;!(l-^2)fc+i 
Combining the above expressions we thus find 



-m'^Nt/{l-fi'-^)_ 



dt 



hk{m) = ^ / (irr2^+3ir^(6r2)e-'^(™'+'''/'"'Vfc(?'V"^) 

Jo 



fe!(l_^2)fc+l 

This form makes it easier to take the microscopic limit. 



dt ^-m'Nt/d-,^) 



(60) 



(61) 



(62) 



(63) 



(64) 



A. Microscopic limit of the bosonic partition function 



The bosonic partition function given by (|24p can now be obtainec 



the Cauchy transform. In the microscopic limit we find (as in 



19|, 



using the expression (j64p for 



25( 1 a factor exp(/i^/2) has been 



removed to ensure independence of ^ for ^ < in the limit fi — > oo) 



+ - 



4/i^/7r 

e^Vm 1 /-i/s^' ds 



duu exp [ 



-m^s— l/4s 



(65) 



^{2mY 2 JO 

The microscopic limit of the Nf = — 1 partition function in the sector of zero topological charge 
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was evaluated in 19|, |25[ with the result, 



X J duueiip[--^\Ko l^-^ j Io{u) + J duuexp[--^\Io yjr^ J Ko{u) . 

= KQ{rn)-e-^^ -^e ^ J duuexp[—^]Koi^]lo{u). (67) 

It is immediately clear that the first term agrees up to an overall normalization constant with the 
first term in the result (j65|) . That also the second terms are in agreement can be seen from the 
substitution of 

i^o(n) = ^^°°f e-*-V4t (68) 
in ()66p which allows the u integration to be performed 

VI. THE PARTIALLY QUENCHED CHIRAL CONDENSATE 

The partially quenched chiral condensate for A^^ flavors is defined by 



1 r„ . / ^Nf,Ni=i 



djln(z';'''"'-\m,,.-- ,mN„m + J\m)] . (69) 
^ ' J=o 



In this section we evaluate the microscopic limit of this expression for Nj = and Nf = 1 and 
derive simplified expressions valid in the strong non-Hermiticity limit. 
For the quenched theory {Nf = 0) we need the ratio 

Z;j^-^'^-^(x|y;/.) _ / fxY^{x^-z]) \ 



J-1 -J I Nf=0 

(y'-x^) (^^J J^N-i{x,y), (70) 



so that the quenched chiral condensate can be expressed as 



VY.Q{m) = 2mHN-i{in,m) , (71) 



m 



where we used the auxiliary kernel defined in (j27p . Below we will evaluate the microscopic limit of 
HN-iim, m). 
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What enters in the chiral condensate for Nf = 1 is the ratio 



Nf=2,N^=l 



(x,mi|y;/i) 



^("ii;/") 



\ y J y y'^-zj 



(72) 



m^pNimi) V y 



— mt 



NN{x,y) pn{x) 
NN{mi,y) PN{mi) 



where the second equahty follows from [171]. The partially quenched condensate for Nf = 1 is the 
condensate as a function of m for fixed physical mass mi. It can be expressed as 



VT,pQ{m,mi) = 2m'HAr(m, m) + 



2m 



m 



mt — m^ 



+ 2m 



PAT (mi) 



1 



m^ — TTIT 



+ 7iN{mi,m) 



■ (73) 



This form makes it easy to take the mi = m limit. The extension to more flavors again only 
requires the function ^{^^{mi^m). 

The integral appearing in the bosonic partition function (|67p has an essential singularity at 



/i = 



25l ]. The function 7i and hence the quenched (17ip and partially quenched (173p condensate 



inherit this non-analyticity. For mi ^ m the partially quenched condensate (see (j73p ) is a function 
of the chemical potential. However, for mi = m it correctly reduces to the derivative with respect 
to m of Zj\f^.=i which is /i-independent in the microscopic limit. 



A. Microscopic limit of the 7i-kernel 

Using the definition of the Cauchy transform, the 'H-kernel can be written as 

N 



HN{mi,m) = y^^Pk{mi)hk{m)/rk 

f 1 ^ 

Jq z^ — m^ ^-^ 



k=0 



To evaluate Ti.]y{mi,m) we substitute 



Pk{mi) = {z -m )[cfc_ipfe_i(z) +Ck-2Pk~2{z) + +pk{zmi/m) 



into dZl]) to get 



1 ^ 

HN{mi,m) = / (fzw{z,z*;fi)^ nyZPk{z^i/'f^)Pk{z*)/rk 

Jc ^ fr^ 



(74) 



(75) 



(76) 
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In the microscopic limit the TL kernel gets multiplied by an additional factor 1/{2N)'^. Using 
expressions for the microscopic limit of the orthogonal polynomials we then obtain 



7i{rni,rh) 



dH 



,2 > 



X / tdt exp( — 2/i t )Iy{ztmi/m)Iy{z*t) . 
/o 



(77) 



B. Strong non-Hermiticity limit 



The expression for 7i simplifies considerably in the strong non-Hermiticity limit, /i ^ 1. Then 
for rhi < m using pi3p we obtain 

|2\ /l-~l2~\ / '• 2 "2 \ 1 

' \z\ mi\ / — m^o\ 1 



ld{rni,rh) = I d z 



Izl'^m'^ 



I. 



exp 



"2 "2 

— m ^2 

1 2 



This is now in a form that can be evaluated analytically. The angular integration gives 



(78) 



7i{rni,rh) = 

The r-integrals are known 
7i{rni,m) 



u-2 



exp 



— rh' 

8/i2 



8/i4m5' 



] r^dr 



r mi 1 Q 

r dr 



(79) 



22( 1 and the final result is 



-1 



— m| 



+ 



m 



9 -^2 ' ^-Q-iy/"-? '-2\ ^-^P 



mj — m 



8/i 



r,2 



77T, 



' 2 

m 



4-1 -TT^ 



(80) 

Inserting this result in (I7ip and (I73p gives respectively the quenched and the partially quenched 
chiral condensate in the strong non-Hermiticity limit. In both cases we need Ti. at equal masses. 
With help of the Wronskian identity this can be expressed as 



2m?Y(m, m) 



u m m 



(81) 



^2/ "2 -2 *2 "2 

In the strong non-Hermiticity limit, the partially quenched chiral condensate also follows from 
the expression of the spectral density as a total derivative. In the quenched case we find 

1 1 



T,Q{m) = I d^z 
Ic 



dz*F{z, z*,'m) + — 
z + rhTT m 



-F{z = m, z* = rh, rh) + 



m 



(82) 



with F given by (see (|136p in Appendix Cp 



19 




FIG. 3: The effect of uiiquenching on the partially quenched chiral condensate is shown by plotting I][/(m) = 
^u.PQi'iTT'i = 0, 2Nra), for one massless flavor with rhi ^ Q, v ~ Q and ji^ ~ 3.84. The strong non-Hermiticity 
result is given by the solid curve. The squares denote the exact microscopic result, and the open circles 
show the result obtained in 26| by matrix diagonalization. 



The subtraction of the singular term amounts to not differentiating ly {fh) except for at m = which 
is compensated for by subtracting the term u/rh. To avoid convergence problems we subtract the 
asymptotic value of the quenched spectral density which is equal to u. The asymptotic contribution 
will be denoted by Sas(?7i). A simple calculation then results in 

2\ / '■ 2 \ /'-2\ /»2 



+ 



m 



4/i2 

rh ' 4/i' 



4/i2 ) (^4/i2 

2 



/ " 2 

/ m 



m 



m m 



le/i^ 



/ ^2 



m 



1 



\ T I 



(84) 



4^2 J ■ ---^ ^4^2 

This result can also be obtained by direct integration the quenched spectral density using polar 
coordinates. The contribution Sas(?7T') is obtained by a direct calculation 



1 



4/i' 



r,2 



dxdy 



1 



X + iy + m 



m 
'4^2 



(85) 



where the j;-integration has been taken over a symmetric interval about x = 0. As it should, the 
result ()84p agrees with (|7ip after inserting (j8ip . 
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Next we study the valence quark mass dependence of the par 



tiahy quenched chiral condensate 



for Nf = 1 and compare with numerical simulations obtained in 26[] for = and nii = 0. In this 
case the result for T,pQ{mi = 0, m) is particularly simple. The chiral condensate is decomposed 
according to (jll2|) 

SpQ(mi = 0, rfi) = T.Q{m) - T,u,pQ{mi = 0, m). (86) 

The quenched contribution was evaluated above (see ()84p ) and is given by (introducing the notation 
u = l/4/i2) 

T,Q{rh) = um + u'^m^[KQ(um'^)lQ{um'^) + Ki{um'^)Ii{um'^)]. (87) 

Notice that this result is only valid when m is inside the domain of eigenvalues. The unquenched 
contribution follows immediately from p26p and (jl27p . or alternatively, can be calculated from the 
strong non-Hermiticity limit of 7i{m, mi) and (1731) . The result for i/ = and rhi = is given by 

tu,PQ{mi = 0,m) = -2ume-'''^'^^Io{m)Ki{um^) + ^. (88) 

The derivation of this result does not require that m is inside the domain of eigenvalues as in the 
quenched case. 

In Fig. [3] we show the result for Tifj^pQ^mi = 0,m) (solid curve) for fi'^ = 3.84. In this 



261 ] by direct 



figure we also give numerical results for an ensemble of 96 x 96 matrices obtained in 
diagonalization of the random matrices. The solid red squares are the results obtained by numerical 
integration of the exact result for the spectral density. 

For complex valence masses we obtain the valence quark mass dependence for v = and nii = 

SpQ(mi = 0,r7i) = T.Q{m) -T,u,PQi'mi = 0,m) (89) 
= 2Mm*e-"™'/2/o(m)Ki(timm*) - 2/m 

+u^mm*rh* {Io{umm*)K(){u7firh*) + Ii{umm*)Ki[umm*)) + um. 

Two three dimensional plots of this result are shown in Figures H] and [5] for fi^ = 3.84 and fP' = 
4 X 3.84 = 15.36, respectively. Note that the x and y axis in the second plot have been scaled 
up by the same factor of 4 as the chemical potential in order that the oscillating area covers 
approximately the same part of the two figures. We observe that the amplitude of the oscillations 
increase strongly with increasing values of jl. 
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Re[E{/,pQ(mi 



1.5, 



-0.5' 



1.1 



Dm 




Ke[m]/2N 



-0.10 



Im[m]/2A^ 



.10 



Q.20 



FIG. 4: The real part of the unquenched contribution to the partiaUy quenched condensate for one massless 
flavor, v — and fi^ — 3.84 as obtained form the strong non-Hcrmiticity approximation. As in Figure 2 the 
X and y axis are scaled by 2N — 96. Note that the result for Im[?fi] = is the same as the full line in Fig. [3) 



It has been shown that the exact microscopic expression for the dependence of the chiral con- 
densate on the dynamical quark mass can be obtained by integrating the unquenched eigenvalue 
density over the complex plane. While the eigenvalue density depends strongly on the chemical 
potential, the chiral condensate does not depend on the chemical potential in the microscopic limit. 
Within a random matrix framework the cancellations leading to this result have been connected to 
the orthogonality properties of orthogonal polynomials in the complex plane and the factorization 
of the /^-dependence of the polynomials in the microscopic limit. The simplifications occur even 
for finite size random matrices and the derivation is simpler in this case than in the microscopic 
limit. The proof of the microscopic limit was obtained by taking the microscopic limit of relations 
derived for finite order polynomials. For completeness, we also have given a derivation of the mass 
dependence of the chiral condensate using complex contour integrations. 

A second effort in this paper was to get a firmer grip on the properties of the spectral density 
that lead to the /i-independence of the chiral condensate. This was achieved in the strong non- 
Hermiticity limit for arbitrary i/, where it was found that an integrability property of the spectral 
density is responsible for the desired cancellations. The interpretation of this property is that 



VII. CONCLUSIONS 
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FIG. 5: To illustrate how the oscillations of the spectral density develop we show (as above) the real part of 
the unquenched contribution to the partially quenched condensate for one massless flavor, v ^ 0, but now 
fi^ — 4x 3.84 — 15.36. The maximal amplitude is about 50 times larger than the scale included in this plot. 

the spectral density has been written as a derivative of the (complex) mass dependence of the 
partially quenched chiral condensate. We have compared the result for the partiahy quenched 
chiral condensate with earlier work, and with exact results obtained by means of complex orthogonal 
polynomials. 

A third result of this paper is the proof of a conjectured identity that was used in the analysis 
of the partition function for one boson. In addition, we have extended this result to arbitrary 
topological charge. 

Finally, we wish to stress that discontinuity of the chiral condensate in the thermodynamic limit 
is due to oscillating terms in the eigenvalue density rather than an accumulation of eigenvalues on 
the imaginary axis: the Banks-Casher formula is not valid for QCD at nonzero chemical potential. 
The nature of the oscillations exemplifies the difficulties encountered by lattice QCD simulations 
at nonzero chemical potential where strong oscillations result in an exponentially suppressed phase 
factor [251] for fi > 

Acknowledgments: We wish to thank Poul Henrik Damgaard for useful discussions. This work 
was supported in part by U.S. DOE Grant No. DE-FG-88ER40388 and by the Carslberg Founda- 
tion (KS). 
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Appendix A. DIRECT DERIVATION OF THE MICROSCOPIC RESULT 

In section IIIII we have manipulated the non-universal finite A'^ expressions into a form which is 
suitable to take the microscopic limit. In this appendix we show that it is also possible to take the 
microscopic limit from the start and deal with the universal expressions throughout the argument. 

The continuum limit of the sum over the index of the orthogonal polynomials is given by 

V F{k) ^ N 2tdtF{Nf). (90) 
fc=o 

The microscopic limit of the orthogonality relation reads 

/ cPzw{z, z*;fi)z* -'Iu{sz*)z-''I^{tz) = 47rfih~^e^^''' d{s - t), (91) 
Jc 

where we used the microscopic limit of the orthogonal polynomials given in (jSSp . The normaliza- 
tion can be verified by means of the large argument asymptotic expansion of the modified Bessel 
functions. 

The microscopic limit of the identity (j2ip is given by 

r ^2 Mz, z*- fi)z* ~-I,{z*t)[z-I,{zt) - m--Umt)] ^ ^ 
Jc z^ -m? 

and the microscopic limit of the expansion (j40p is given by 

z-''l{z)-m-''h{m) 



_g2 _ ^2 



[ A{rh,s)z~''I^{zs)sds, (93) 
Jo 



with expansion coefficients equal to 

A{m,s) = Ky{ms)Iu{ra) — Ky[m)Iy{ms). (94) 

Using the definition of a derivative we obtain 

m^'^[I' {m) /^(m)] = 2m / A{m, s)m~'^ Iu{rhs)ds. (95) 

m Jo 

The microscopic spectral density for Nf = 1 is given by P, Q] 

(u) . w(^^*]fj) ^^^2^2\Iu{z*t)U{zt) Iy{zH)I^{rht) U{z) 







\z\'^'^ lui'lTT-) 



(96) 



By adding and subtracting Iu{mt) to Iu{zt) and luirh) to Iu{z) the spectral density can be rewritten 



as 



z*" 

Iy{rht) 



(y) . wyz,z , 

p^^^{z,z ,m;/i) = 



X [z-^Uzt) - m-^h{mt)] - [z-^h{z) - m-/,(m)]^^ . (97) 

Iu[m) J 
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The chiral condensate is given by 

S(?n) = / (fz p^^l_^{z,z*,m;fi). (98) 

Because of the expansion ()93p and the orthogonaUty relation (|9ip . the first bracketed term in Eq. 
([97]) does not contribute to the chiral condensate resulting in 

S(m) = rn[ d^z^^^M^ dtfe-^A^^^ " ^7^-^^^ ^^^'^ y,^',^ ■ (99) 

Jc 27r/i Jo z — m z*^ Iv{iti) 

Finally, after inserting the expansion (|93p and applying the orthogonality relation (j9ip . the deriva- 
tive relation (j95p results in 

lu[m) m 

which is the correct /^-independent chiral condensate after including the contribution from the zero 
modes. 



Appendix B. CHIRAL CONDENSATE FROM COMPLEX CONTOUR INTEGRATIONS 

In this appendix we derive the chiral condensate by integration over the eigenvalue density using 
complex contour integrations. We only work out the case of = 0. 
Let us start by writing the chiral condensate for Nf = 1 as 

±{m) = [ dt±{rh,t) (101) 

where 







t{m,t) = [ (i^z— 3— r^^^_^(i,z*,t;77i,/i), (102) 
Jc z + m "f ' 

and 

pg;._,(z, z*,t- rh, A) = -^e-(^^+^")/«A-K„(|l| V4A^)te-^^^*^/o(r t)/o(mO " ,(103) 



which follows from (j96p after use of (j92p . Since the chiral condensate is real, there is no need to 
consider imaginary contributions to T,{m,t). 

We write z = x + iy and do the integral over y by a contour integration. Since we perform the 
integration over y by a contour integration, we have to distinguish the cases x < —rh and x > —rh 
and decompose Io{z) as 

I^{z) = -^{Ko{z)-Ko{-z)). (104) 
vr J—z 
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For /o(i*t), a similar decomposition is used. The pole is at ?/ = i{x + rh) so that 



-m and I* — > 2a; + m. 



(105) 



Completing the contour integral we obtain for the pole contribution 



S(m, t) 



—rht 



dx{2x + rh)e 



2i-t2H§i±2H_ 2^2*2 m(2x + m) lo(mt) 



X [e{x + m)[-isgn(2x + m)i^o((2x + m)t)Io{m) + iKo(m)/o((2x + m)t)] 

-d{-x - 7n)[isgn(22; + m)Ko(-(2x + m)t)Io{m) - iKo{-m)Io{{2x + m)t)]] (106) 

There is also a contribution from the jump across the cuts along the imaginary axis. However, 
since these contributions are purely imaginary we do not further analyze them. 

To disentangle the different contributions, we rewrite the hyperbolic Bessel functions as 



Koix) = Koi\x\)-7TiIoi\x\). 



(107) 



After the cancellation of two terms for x + ?n < we observe that the integrand is given by the 
same expression as for x + m > 0. Using y = 2x + ?fi as new integration variable the expression 
for Re(S(m,t)) simplifies to upon integration, 



Re(S(m,t)) 



^ / dyye ^ 



0{y)h 



m\y\ 



4/i 



r,2 



(/o(m) sgn(y)i^o(|y|0 - KQ{m)lMt)) 



+0{-y)Ko 



m\y\ 



Io{rh) sgn{y)Ioi\y\t) 



(108) 



The integrals can be simply rewritten as 
Re(S(m,t)) 



mt Io{mt) -2A2t2-|4 ^ -j4 
^ / dyye 

Jo 



2/i2 /o(m) 



(109) 



[lo Koiyt) +Ko(^] Io{yt)Mm) - lo ( ^ ] Io{yt)Ko{m) 



my 
4fj^ 



my 
4fj^ 



The integral of the terms multiplying /o(fTi) was evaluated in 



the last term in the above equation can be found in 

Io{rht) 



25( 1 whereas the integral containing 



2M- We finally find 



Re(E(m,t)) = 2mt 



hijh) 



KQ{mt)lQ{m) - Ko{m)Io{mt)] 



(110) 



in agreement with results derived in subsection IIII B 11 
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Appendix C. SPECTRAL DENSITY AS A TOTAL DERIVATIVE 

In this appendix we obtain the chiral condensate by writing the spectral density as a total 
derivative. This will be worked out in the strong non-hermiticity limit only, first for i' = ^ and 
then for general z^. 



A. Strong Non-Hermiticity 

The microscopic limit of the spectral density (j29p for Nj = 1 is 5, ol] 

p^;;]^,{z,z*,rn;f^) = J-L^, ^LI_j (111) 
x( [\tt e-^f"'"' I^{zt)I^{z*t) - C dt t Q-"^^^'^ I^{rht)I^{z*t) 

The term with the first integral is the quenched density and the second term gives the effect of 
unquenching, which motivates the notation 

It is easily checked that pNf=i vanishes at z = m since pQ = lim^_^£ pjj. This trivial observation 
which holds since the fermion determinant vanishes when an eigenvalue is equal to the quark mass 
will be very useful below. 

In the limit of strong non-Hermiticity p? ^ 1 with x/{2p?') < 1 and (x + rfi)/(4/i^) < 1 such 
that the saddle points of the i-integrations are inside the interval [0, 1], the integrals in (|111|) can 
be approximated by 

-1 />oo 

dtte~^'^^^^Iu{ix-iy)t)Iu{mt) ^ / dtte~'^^^^^ Iy{{x - iy)t)lQ{rht) 
Jo 

^2 j_ vj,2 



1 +m / fj^ix - ill) , , , 



Within this limit the quenched spectral density for topological charge u reduces to 

j^) ^ '^^2*K^{uzz*)h{uzz*), (114) 

while the unquenched part reads 

>) = '^2z*K,{uzz*)U{umt)^4^^e<'^"-'^"^l^ (115) 
vr lv{rn) 

where we introduced the abbreviation u = 1/4/i^. Below we will express p as a total derivative 
and show that the resulting chiral condensate is independent of p. Because the chiral condensate 
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is independent of fi we recover the full microscopic chiral condensate even though we work in the 
limit of strong non-Hermiticity. Before we present the argument for general topological index we 
first work through the simpler case u = 1/2. 

1. The case of v ~ ^ 

The microscopic spectral density for = ^ and Nf = 1 is given by 

P^N^lliz.z* ,m;fi) = pQ~^\z,z*;fi) - p''^' ^\z, z* ,m; ft) (116) 
with (using the notation u = 1/4/i^) 

^(|) _ ^^-^i^--uii* + ^m^ e —e ,urhz* _ ^-urhz*\ (117) 



and 



p[|)= lim_pc; = -(l-e-2--*). (118) 
It is straightforward to integrate pjj with respect to z*: 

pif = -9,4 ^^ ' ^ — , i] + Ld,4 ^^ '^'"^^ — , —] (119) 

vr m — z TT z + m 

with 

The limit rh ^ z and the derivative with respect to z* can be interchanged for smooth functions. 
However, taking the limit first can generate singularities that result in spurious contributions after 
differentiating with respect to z* . Taking the limit m — > z in (|119p by naive differentiation leads 
to 



cosh z e _|_ \ 

uz + uz r-r-r + 



TT [ sinhz 2z 

To eliminate the spurious singular terms we subtract 



(121) 



pf = --dA^--\ (122) 

^ vr smhz z 

from the quenched spectral density. The chiral condensate can then be expressed as 

S(m) = j dxdy^^ipf - pf - pf). (123) 
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The contribution of the unsubtracted derivatives vanishes at z = rh. After partial integration we 
thus find 

^ coshm 1 h/2{m) 

S m) = —-: - = - — — . (124 

smhjTT, m -'i/2("^) 

After adding the contribution from the zero modes given by ly/m we find 

h/2{m) 1 
him) = — - H 

= (125) 
in agreement with the result for arbitrary topological charge to be discussed next. 



2. The general v case 



We now turn to the case with arbitrary u. Similar to the case = \ above, we express p[J as 

. (126) 

with 

F{z,z*,m) = 2ze'''^'^^~^^^'^^^[mz*uK^{urz*)ll{umz*) - zz*uI^{umz*)Kl{uzz*)] - 2z, 

Iy[m) 

= FAiz,m)G{mz*,zz*) -2z. (127) 

Here, 

FA{z,m) = 2iexp(n(m2 - z'^)/2)I^{z)/I,y{m), (128) 

and 

G{mz* , zz*) = mz*uKy{uzz*)l'j^[urhz*) — zz*uIp{umz*)Kl{uzz*). (129) 

The term 2z has been subtracted to eliminate singularities at z = itm. Using the Wronskian 
identity 

ll{x)K,{x) - Iu{x)K{x) = i (130) 

one easily shows that 

G{mz\zz*)FA{z,m)\^^^, = 2z. (131) 
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Therefore, 



so that 



Fa{z, 'm)G{mz* , zz*) — 2z = {w? — 2p')h{rh, z, z*) 



(132) 



d^[FAG-2z] 



z=z''=m 



-d,[FAG-2z 



z=z''=m 



(133) 



with h{m, z, z*) a function that is regular at m = z. 
The quenched part of the spectral density is given by 

Pq\z,z*;i1)= ]:iT[ipp{z,z*,m]ji). 



Naively interchanging this limit with dz* results in 

Ml- -*\ _ 1 (a FA{z,z)drhG{mz,zz* 



Pq i^^z ) = - 



dz 



G{mz, zz*)drhFA{z,m)\^. 



(134) 



(135) 



2z 2z 
However, the second term may give rise to singular contributions that would be absent if we would 
have differentiated with respect to z* before taking the limit. Therefore these contributions have 
to be subtracted so that the quenched spectral density is given by 



H.. _ 1 drh[FA{z,rh)G{mz, zz*)]j 



IT \ 2z 
The chiral condensate can be written as 



Sf.Sing 



G{rhz, zz*)drriFA{z, m) 



2z 



cfz 



1 1 



m + Z TT 



drh[FA{z,m)G{7hz, ZZ*))]r 



2z 



(136) 



(137) 



-Sing 



G{mz, zz*)dmFA{z,m)\^. 
2z 



FAiz, z)G{rhz, zz*) - 2z 



rh?' — 



After partial integration with respect to z* we find 

'drh[FA{z, z)G{mZ, ZZ*)]ra=z 



Il(m) 



2z 



(138) 



-Sing 



G{fnz, zz*)drhFA{z,m)\^, 
2z 



+ 



d^[FA{z,z)G{mz,zz*) - 2z] 
2m 



z=z'=—m 



After using (jl33p only the singular terms remain. Because of (jl32p there is no singularity at z = 
after using (jl33p 



T,(rh) = Sins 



G{rhz, zz*)drh Fa{z, m) 
2i 



2=2*=— m 



ly {m) rh ' 



(139) 



which gives the correct result after including the contribution from the zero modes. 



[1] P. E. Gibbs, Glasgow University preprint PRINT-86-0389 GLASGOW (1986), unpublished. 



30 



[2] D. Toublan and J. J. M. Verbaarschot, Int. J. Mod. Phys. B 15, 1404 (2001) |hep-th/0001110| . 

[3] I. Barbour, N. E. Behilil, E. Dagotto, F. Karsch, A. Moreo, M. Stone and H. W. Wyld, Nucl. Phys. B 

275, 296 (1986). 
[4] T. Banks and A. Cashcr, Nucl. Phys. B169, 103 (1980). 
[5] J. C. Osborn, Phys. Rev. Lett. 93, 222001 (2004) iarXiv:hep-th/0403131| . 

[6] G. Akemann, J. C. Osborn, K. Sphttorff and J. J. M. Verbaarschot, NucL Phys. B 712, 287 (2005) 



[arXiv:hep-th/0411030| . 

[7] J. C. Osborn, K. Sphttorff and J. J. M. Verbaarschot, Phys. Rev. Lett. 94, 202001 (2005). 

[8] L. RavagU and J. J. M. Verbaarschot, Phys. Rev. D 76, 054506 (2007) farXiv:0704.1111l [hep-th]]. 

[9] E. V. Shuryak and J. J. M. Verbaarschot, NucL Phys. A 560, 306 (1993) |arXiv:hep-th/9212088| . 
[10] J. C. Osborn, K. Sphttorff and J. J. M. Verbaarschot, Phase Diagram of the Dirac Spectrum at Nonzero 

Chemical Potential, in preparation. 
[11] J. J. M. Verbaarschot, Phys. Rev. Lett. 72, 2531 (1994) [arXiv:hep-th/9401059| . 
[12] J. J. M. Verbaarschot and T. Wettig, Ann. Rev. Nucl. Part. Sci. 50, 343 (2000) |hep-ph/0003017| . 
[13] K. Sphttorff and J. J. M. Verbaarschot, Nucl. Phys. B 683, 467 (2004) jarXiv:hep-th/0310271| . 
[14] G. Akemann, Acta Phys. Polon. B 34 (2003) 4653 jhep^th/0~3'07lT6] . 
[15] G. Akemann and G. Vernizzi, Nucl. Phys. B 660 (2003) 532 [hep-th/0212051| . 
[16] M.C. Bergere, [hel>th/0311227| . 
[17] M.C. Bergere, [hep-th/0404126l . 

[18] G. Akemann and A. Pottier, J. Phys. A 37, L453 (2004) 'math-ph/0404068| . 
[19] K. Sphttorff and J. J. M. Verbaarschot, Nucl. Phys. B 757, 259 (2006) [arXiv:hep-th/0605143| . 
[20] J. C. Osborn, PoS LAT2006, 142 (2006) |arXiv:hep-lat/0610016| . 
[21] K. Sphttorff, J. J. M. Verbaarschot and M. R. Zirnbauer, |arXiv:0802:2660l [hep-th]]. 
[22] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products (sixth edition) Academic 
Press (2000). 

[23] V.L. Girko, Theory of random determinants (Kluwer Academic Pubhshers, Dordrecht, 1990). 
[24] M.A. Stephanov, Phys. Rev. Lett. 76 (1996) 4472 [arXiv:hep-lat /9604003| . 

[25] K. Sphttorff and J. J. M. Verbaarschot, Phys. Rev. D 75, 116003 (2007) |arXiv:hep-lat /070201 1| . 
[26] M. A. Halasz, A. D. Jackson and J. J. M. Verbaarschot, Phys. Rev. D 56, 5140 (1997) 
|arXiv:hep-lat/9703006| . 



